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This paper explores features of an idealized mathematical machine (algorithm) that would be 
capable of reconstructing the gravitational nature (the multipolar structure or spacetime metric) of 
a compact object, by observing gravitational radiation emitted by a small object that orbits and 
spirals into it. An outline is given of the mathematical developments that must be carried out in 
order to construct such a machine. 
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I. INTRODUCTION 

Gravitational wave detectors such as LIGO (the 
Laser Interferometer Gravitational Wave Observatory) 
are rapidly increasing their sensitivity, making precise 
measurement of waveforms emanating from massive grav- 
itating objects a reality in the near future. The launch of 
LISA (Laser Interferometer Space Antenna) will further 
increase our observational capacity 

One method of mapping out the space time of strong 
field regions is observing the waveforms of extreme and 
intermediate mass ratio inspirals (EMRI's and IMRI's). 
The physical scenario is presented in Figure [T] A low- 
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FIG. 1: Caricature of an EMRI 

mass inspiraling object (a "probe", e.g., a particle or 
carrot or neutron star ) moves through the background 
spacetime on nearly geodesic orbits around a central com- 
pact object. As it does so it samples the geodesic struc- 
ture of the background manifold, as warped by the com- 
pact object, and broadcasts this information via gravita- 
tional radiation to detectors such as LIGO. The intrin- 
sic periods associated with the particle motion, and the 
change of these periods as the particle inspirals give us a 
way of characterizing the spacetime in which the particle 
is moving. 

The current mathematical formulation of this problem, 
waveform generation and data analysis techniques, is suf- 
ficient only for "observing" Kerr EMRI inspirals 0, 0, H| ■ 



In effect, current techniques presuppose that the central 
object is a Kerr black hole, i.e. the unique Kerr space- 
time favored by the no-hair theorems @. In other 
words, it is assumed that the axioms of cosmic censor- 
ship and causality hold and that the only parameters to 
be determined are the mass and spin of the black hole 
and orbital parameters of the probe. 

Suppose, however, one would like to entertain the idea 
that something more exotic @, H, Q may be created in 
the strong field regions of the universe, and would like 
to find a way of observing, rather than presupposing, 
what these features are. Suppose that one would like to 
put ideas such as cosmic censorship and causality to an 
experimental test rather than using them to aid the data 
analysis. How would one go about, in effect, drawing 
the bunny out of the hat in Figure [TJ by watching the 
radiation emitted by the inspiraling object? 

The outline of a mathematical machine, although com- 
plex, that could possibly do so, is the subject of this pa- 
per. Subsequent papers in this series, Eqi, m 0, win 

develop tools that may make possible such a machine. 

In Sec. |TT] we comment on the existing concept for 
EMRI searches in gravitational wave (GW) detector data 
and highlight at each step the mathematical features that 
make it possible. 

In Sec. IHII we give a suggested formulation of the prob- 
lem that could, in principle, encompass, as the central 
object, all stationary axisymmetric vacuum spacetimes, 
and we highlight ideas from integrable systems and exact 
solutions of the Einstein equations that could underpin 
the desired machine. 

Finally, in Sec. IIIII we also identify crucial points in 
the mathematical understanding of the field equations 
that must be sorted out in order to make such a machine 
viable. 



II. THE EXISTING EMRI WAVEFORM 
GENERATION MACHINE 

Current calculations of EMRI waveforms have been re- 
stricted mainly to inspirals around Kerr black holes. 

A schematic sketch of the current waveform generation 
technique, as implemented by Drasco and Hughes 0], and 
a search algorithm are given in Figure f2] 
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FIG. 2: EMRI waveform generation machine 



In step A it is assumed that the probe particle is mov- 
ing around a Kerr black hole whose mass and spin are 
known. The motion of the probe particle results in a per- 
turbation on the background spacetime which is trans- 
lated using Teukolsky's equation 0, [l3| into a waveform 
in the asymptotic region where the detectors are located, 
step B. The observational step C involves detection, and 
estimation of the source's parameters, via matched filter- 
ing. The noisy nature of the data of all GW detectors 
makes this final step a necessity: a method for differen- 
tiating between features of the detected signal that have 
their origin in the EMRI signal and those that have their 
origin in detector noise is required. Once a detection is 
made and the parameter estimation step for the mass, 
spin and orbital parameters of the inspiraling object is 
conducted, we attribute the gravitational wave event to 
having been triggered by an inspiral around a Kerr black 
hole. The uniqueness of this identification, is in most 
cases, assumed. 

Let us now analyze each step of the process in greater 
detail and note the features of the calculation that make 
it tractable. The Kerr spacetime is of Petrov type D (for 
the definition of Petrov type see [T3|) and admits a full set 
of isolating integrals, or constants of motion; namely, the 
rest mass /x of the particle, energy E, axial angular mo- 
mentum L z , and the Carter constant Q. While the first 
three constants can be trivially obtained from the met- 
ric symmetries, Eq. (TJJ), the Carter constant Q is more 
subtle. Discovered by separation of the Hamilton Jacobi 
equation by Carter in 1968 [15| , it plays an important 
role in every step of the calculation. In step A the con- 
stants of motion (or action variables, as they are known 
in the field of dynamical systems) give us the power to 
fully describe the orbit: the action variables uniquely 
identify the orbit of the test particle around the compact 
object and describe its physical confinement. The angle 
variables identify where on the trajectory the particle is 
located at a given time, while traversing the orbit. 

In step B the particle's motion in its orbit serves as 



a source of the Teukolsky equation, which is used to 
translate the perturbations caused by the particle in the 
strong field region into the gravitational waves we ob- 
serve with our detectors. The Teukolsky equation can 
be solved by means of separation of variables (HI; the 
analysis performed by Teukolsky is only valid in Petrov 
type D spacetimes. It is this separability feature of 
the Teukolsky equation that is exploited by Drasco and 
Hughes 4] to perform the translation of particle motion 
from the strong field region into gravitational waves in 
the asymptotic region. There turns out to be a very 
deep relationship between the separability of differential 
operators such as those governing the Teukolsky equa- 
tion, and the existence of a second-order Killing tensor 
on Petrov type D spacetimes [l6| (Ch 5 and the refer- 
ences therein). In the Kerr spacetime, the Carter con- 
stant can be attributed to the existence of second-order 
Killing tensor, and the Teukolsky equation's separability 
in Boyer-Lindquist coordinates can be seen as a natural 
result. 

In step C the waveforms computed for Kerr in steps A 
and B are used as templates for matched filtering. The 
effect of the Carter constant is evident in these templates, 
in that the Fourier spectra of the waveforms themselves 
are built up of harmonics of three fundamental frequen- 
cies that characterize the orbit (the three frequencies map 
directly onto E, L z , Q). The matched-filtering step is 
essential in noisy experimental environments. However, 
matched-filtering is also very limiting, in that it requires 
that we know the form of the templates before performing 
data analysis. As a result, the current EMRI template 
bank consists mainly of Kerr inspirals. 

Other isolated examples of explorations of spacetime 
mapping, or the reconstruction of the multipole-moment 
structure of the central object from observed gravita- 
tional radiation, include Ryan's exploration into the fea- 
sibility of detection of the multipolc moments of a boson 
star [8j. Collins and Hughes |l7[ provide an important 
contribution toward the formalism for mapping space- 
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time around "bumpy blackholes" by studying the orbits 
confined to the equatorial plane around "bumpy" static 
objects described by the Weyl class of solution to the 
vacuum Einstein equations. 

The matched-filtering approach hampers observation 
by requiring that one initially compute the model of the 
central object, work out possible inspiral templates from 
the resulting spacetime, and then conduct the match- 
ing and parameter estimation. An effective spacetime 
mapping algorithm, beyond Kerr spacetimes, would re- 
quire an enormous number of templates and even then 
one could not possibly hope to cover all possible scenar- 
ios. In effect, the observational power of the LIGO and 
LISA detectors in the EMRI inspiral problem is limited 
by the models we can conceive and calculate. To date, no 
general framework exists that will allow us to effectively 
map the strong field region around an unknown object. 

In the next section we formulate the EMRI problem 
for general stationary axisymmetric vacuum spacetimes, 
and we suggest methods in which the ideas of integrable 
systems can be applied to make a general detection algo- 
rithm possible. 

III. FORMULATION OF THE EMRI PROBLEM 
FOR AXISYMMETRIC STATIONARY 
SPACETIMES 

Consider a central body with arbitrary multipole mo- 
ments and a probe particle moving in the vacuum space- 
time around it. 

The line element of this spacetime can be represented 
in the Lewis Papapetrou form: 

ds 2 = e~ 2 ^ [e 2 ~<(dp 2 + dz 2 ) + R 2 dtf] - e 2 ^(dt - tud^ 2 , 

(1) 

and is entirely determined by solutions of the Ernst equa- 
tion [Lsj for the complex potential £, 

»(£) V 2 £ = V£ ■ V£ (2) 

Any axisymmetric, stationary vacuum solution can be 
identified by means of a bi-infinite sequence of numbers 
physically interpreted as multipole moments Mi and Si. 
This corresponds to choosing a particular element of the 
Geroch group. As first conjectured by Geroch 0, [2(| 
and later proved b y H oenselaers . Kinnersley, Xanthopou- 
los and by Chitre [U, [H [H HI H| in the late 1970's in 
a series of papers that lead to the HKX transformations, 
these numbers uniquely identify a spacetime. Once they 
are known, the spacetime is in principle determined. Sub- 
sequently, a number of other solution-generating tech- 
niques have been developed that allow one to determine 
the explicit form of the functions in the metric |T]) by 
mapping a given solution of ([2]) onto another [26| . 

We speculate that it may be possible to exploit these 
mappings to help develop an algorithm that could in prin- 
ciple limit, if not completely determine, the multipole 



structure of a spacetime from its EMRI inspiral waves. 
We further speculate that it may be possible to do so 
without resorting to matched filtering and its need for a 
priori guessing the structure of the central object. 

There are several uncertainties implicit in expanding 
the model shown in Figure [5] from Kerr to a general 
method for mapping spacetime, all of which have to be 
addressed before a spacetime reconstruction algorithm 
becomes practical and before we can determine how much 
information can, in practice, be gleaned from an EMRI 
inspiral event. These uncertainties include: A) whether 
or not an explicit action-angle variable prescription can 
be found that gives us access to the full description of 
the probing particle's geodesic orbit. B) If A) is indeed 
possible in large regions of more general spacetimes, is 
it feasible to attempt to explore the perturbation prob- 
lem on a general background and what form would that 
calculation take. C) Observationally, one is only privy 
to partial knowledge of the gravitational wave emission 
originating from a sequence of geodesic "snapshot" or- 
bits (portion of the orbit in which the radiation-reaction- 
induced evolution of integrals of the motion is negligible), 
as in Part B of Figure [5] The noise of the detector, the 
effect of the mass of the probing particle on its motion 
through the background spacetime that moves it off the 
geodesic trajectory (self force) and the length of the ob- 
servation or validity of the adiabatic approximation all 
conspire to complicate the signal. What is sought is a 
method to extract the signal from the noise and a sub- 
sequent representation of that signal that allows one to 
clarify the nature of the non-geodesic effects and quan- 
tify which parameters describing the central object can 
be obtained with certainty. All of these questions and 
uncertainties are addressed in the next three paragraphs. 



A. Orbital description 

In many calculations in General Relativity it is implic- 
itly assumed that it is possible to find four constants of 
geodesic motion that describe the orbital motion of a test 
particle, and it is the bias of the author that it is indeed 
so. 

However, in the past, such intuition has proven faulty. 
Poincare's study of the three-body problem and the ad- 
vent of our understanding of deterministic chaos forever 
banished the ideal of finding a beautifully simple closed 
form description of particle motion in Newtonian grav- 
ity (27[. The Henon-Heiles problem warns that, in the 
event of the system being chaotic, perturbation theory, 
while it will yield a computational result, will fail to ac- 
curately represent the dynamics of the system [28|]. Nu- 
merical studies into the orbital nature of SAV spacetimes 
are conducted i n 12011 and are discussed further in Paper 
II of this series [10( ■ Understanding the interrelationship 
between the existence of Killing tensors and the integra- 
bility properties of the Ernst equation and solution gen- 
eration techniques may provide a constructive method of 
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finding the invariants in question; see Papers III and IV 
in this series [111, [l2j . Such an investigation may also shed 
light on a possible approach to the resolution of problem 
B) 



B. Translation of orbital motion into detection 
region. 

The translation of the effect of particle motion in the 
strong field region into waves in the asymptotic region 
and subsequent coupling to the detector requires that 
one solves the perturbation problem off of all SAV so- 
lutions. Although this task may seem daunting, and 
very little appears to have been done on perturbations 
off more general SAV backgrounds than Kerr, the idea of 
GW travelling outward can, in some sense, be viewed as 
a particle perturbation traveling along a series of plucked 
strings or geodesies toward the asymptotic region. Just 
as in the Petrov type D case, the integrability properties 
of the wave equations (or perturbation equations) in the 
general spacetimes should be related to the background 
geodesic structure of the spacetime one is perturbing off. 
The two aspects of the Teukolsky analysis [13^ that make 
the problem tractable in type D spacetimes, namely the 
decoupling of the perturbation equations and separabil- 
ity have subsequently been understood in terms of second 
order Killing tensors [l6| ■ The extension of this work to 
higher order Killing tensors would provide a point of en- 
try to solving Problem B). 



C. Detection and Noise Control 

A possible scenario in which the matched filtering cri- 
terion could be lifted is by finding an experimental re- 
alization of the solution method for the Ernst equation 
employed by [3(| [M| initially given in [32j . What is done 
is to introduce a linear potential matrix $ much like a 
wave function in quantum mechanics. In the equation 
governing $ the Ernst potential enters as an unknown 
field, a quantum mechanical potential well per analogy. 
A great deal is known about the properties of <3? [30l | . 
without a priori specifying the explicit gravitational field 
or Ernst potential involved. In effect, the <!> serves as a 
carrier or equivalence class for the gravitational poten- 
tial being observed allowing its properties to be known 
without specifying the entity itself. 

An example of using the known properties of the so- 
lutions of an equation to aid detection within a noisy 
environment, without explicitly modeling the waveform 
and thus knowing the source of the waves, can be found 
in the form of the KdV equation describing shallow water 
waves |33j. This example is much simpler than the SAV 
problem and is accompanied by tangible physical inter- 
pretation. Furthermore, the KdV equation shares many 
of the same mathematical properties as the GW, SAV 
problem, and I speculate that it may be worthwhile ex- 



ploring to provide insight on how to proceed in the GW 
case. The KdV equation admits solution by means of 
the inverse scattering method described in the previous 
paragraph, albeit much less complicated [3(| [Hj]. One of 
the features of the solution identified in the analysis is 
the dispersion relation. If one is unfortunate enough to 
attempt to detect shallow water waves on a pond while 
one's child is splashing in the foreground of the machine, 
filtering the data using this dispersion relationship, and 
some knowledge of the functional form of the waves origi- 
nating further out may suffice in removing the child from 
the measurement without ruining a good day's play. How 
exactly, to effect such a filter for a gravitational wave ex- 
periment is at present unclear. Two things, however, 
are certain: we cannot remove the noisy child and, with- 
out the filter, our observational power in the context of 
spacetime mapping is limited. 

For our SAV gravitational-wave problem, an approach 
to the representation of the GW data could be as follows. 
Observe that the solution to the geodesies equations of 
Kerr can be written down in terms of Weierstrass's ellip- 
tic functions |34j, the poles of which are related to the 



constants of motion, or Killing tensors. Identify segments 
of the real waveform with the poles associated with the 
corresponding snapshot waveform. This will allow us to 
compute pole tracks as the particle moves from geodesic 
to geodesic. The self force calculations should provide 
the theoretically expected track. In a general spacetime 
mapping setting, make a similar identification and sup- 
pose, initially, that the inspiral is around a Kerr object. 
Should the pole track begin to deviate from Kerr, sys- 
tematically adjust the lower order multipole moments of 
the model so that the observed track best matches the 
theoretical one. For this the self force calculation on a 
general background would be required. 

An expression of the metric that is explicitly 
parametrized by means of the multipole moments will 
facilitate calculations. In practice, for static spacetimes 
which are a member of the Weyl class [35] , a metric which 
is parametrized by means of all the mass multipoles al- 
ready exists. For SAV spacetimes an explicit form of 
the metric is not known, but known special cases may 
be helpful in developing our proposed techniques: The 
Manko-Novikov spacetime [361 ] provides an explicit met- 
ric in which all mass moments and some, but not all, 
the spin moments appear explicitly in the metric. Other 
metrics that would be of astrophysical interest and for 
which the Ernst potential is explicitly known are those 
of a compact object surrounded by a disk (30j . 



IV. CONCLUSION 

There will always remain experimental uncertainty as 
to how well one can determine the structure of the central 
object and thus the extent to which one can verify the 
validity of the no-hair theorems || or confirm the exis- 
tence of more complex central objects. Ironically, quan- 
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turn mechanics and the act of measurement itself force 
us to play dice in determining the details of Einstein's 
theory in practice. 

Possibly the experimental and data analysis challenge 
is to find the means by which we can learn the most. 
If the relationship between curvature content of the 
Weyl Tensor, as encoded in the Ernst potential, and the 
geodesic structure, can be understood in detail and ap- 
propriately exploited, it may lead to a powerful experi- 
mental application of the mathematical development in 
the field of exact solutions over the last few decades. A 
possible framework by which this can be done has been 
suggested in this paper. 

This paper is presented as a question about the feasi- 
bility of this mathematical machine. You are encouraged 



to find its flaws. In subsequent papers [l(| [0, EH > some 
of the ideas presented here will be placed on a firmer 
mathematical footing, thus laying the foundation for the 
construction of this machine. 
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